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We find a new type of optical vector soliton that originates from trapping of a dipole mode by a 
soliton-induced waveguide. These solitons, which appear as a consequence of the vector nature of 
the two component system, are more stable than the previously found optical vortex-mode solitons 
and represent a new type of extremely robust nonlinear vector structure. 
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Perhaps one of the most desirable goals of Optics is the 
development of purely optical devices in which light can 
be used to guide and manipulate light itself. This moti- 
vation explains the growing interest in self-guided beams 
(or spatial optical solitons) and the recent theoretical and 
experimental study of spatial solitons and their interac- 
tions (l). It is not only the reconfigurable and steerable 
soliton-induced waveguides created in a bulk medium 
that are of particular practical interest, but also a spatial 
soliton that guides another beam (of a different polariza- 
tion or frequency). This may become a completely new 
object, a vector soliton, with an internal structure and 
new dynamical properties which yield surprising results 
for the stability of such an object even in simplest cases. 

Complex phenomena induced by the vector nature of 
nonlinear wave equations arise in many fields of Physics 
and are already firmly placed in the realm of condensed 
matter physics, the dynamics of biomoleculcs, and non- 
linear optics. Recently an interest in these structures 
and the theoretical possibilities they offer has been re- 
newed because of their experimental realization in differ- 
ent physical contexts. For instance, vector phenomena 
have been observed in Bose-Einstein condensation, with 
vortices in multicomponent condensates J^] or in nontriv- 
ial topological defects due to interspecies interaction || . 
Finally, vectorially active optical media are also being 
currently investigated because of many new characteris- 
tics they provide, as compared to scalar systems j|. 

In this Letter we study a vector object formed by two 
optical beams that interact incoherently in a bulk non- 
linear medium. If the nonlinearity of the medium is self- 
focusing, an isolated beam, under proper conditions, will 
form a self-trapped state - a spatial optical soliton |Q. 
Such a soliton changes the refractive index of the medium 
and creates a stationary effective waveguide. A second 
beam of a much lower intensity is subjected to the in- 
duced change of the refractive index and can be trapped 
as a localized mode of that waveguide. 

From linear optical waveguide theory, we expect that 
a radially symmetric waveguide can host different types 
of modes with more elaborate geometries [Figs, l(a-c)]. 
However, at higher intensities of the trapped beam one 



must regard the two beams as components of a vector 
soliton, self-trapped by a self-consistent change of the 
refractive index induced by both beams. In this case we 
cannot treat the shapes of the beams as independent and 
it is not trivial to conclude whether we may obtain states 
which are a composition of a lowest-order state in one 
component and a high-order state in the other one. 
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FIG. 1. Intensity distribution of (a) fundamental, (b) vor- 
tex, and (c) dipole modes supported by a radially symmetric 
fundamental soliton-induced waveguide. 

Recalling a previous work on the existence and stabil- 
ity of two-component one-dimensional solitons ||, one 
may expect to find at least two types of such complex 
objects in the two-dimensional case. The first type is 
the recently discussed two-dimensional vector soliton || 
which has a node-less shape [e.g., as shown in Fig. 1(a)] 
in the first component and a vortex on the second one 
[Fig. 1(b)]. The second type introduced here maintains 
a similar shape for the first component, while the sec- 
ond beam develops a node along a certain direction [Fig. 
1(c)] forming what we call a dipole-mode vector soliton. 

The purpose of this Letter is two- fold. First, we dis- 
cuss the stability of the vortex-mode vector solitons and 
show that these objects are linearly unstable and can de- 
cay into dipole-mode vector solitons. Secondly, we prove 
that these dipole-mode solitons exist for an ample range 
of relative intensities of their components and show that 
they survive both small and large amplitude perturba- 
tions, their propagation dynamics resembling that of two 
spiraling beams |7j . 

We would like to emphasize here that both results are 
highly nontrivial. While it is commonly believed that 
asymmetric solitary waves possess a higher energy, and 
should be a priori unstable, our results demonstrate that 
the opposite is true: An excited state with an elaborate 
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geometry may indeed be more stable than a radially sym- 
metric one and, as such, would be a better candidate for 
experimental realization. We stress here that the recently 
discovered method of creating multi-component spatial 
optical solitons in a photorefractive medium || would 
allow a simple and direct verification of our theory, in- 
cluding the questions of soliton existence and stability. 

The outline of the Letter is as follows. First we for- 
mulate the model for our system. Next, we proceed with 
the study of vortex- mode vector solitons in the reasonable 
parameter range of the model. We recall previous studies 
on the issue of existence of such states |^|| and add new 
results arising from a linear stability analysis and numer- 
ical simulations of the dynamics of these solitons. Having 
concluded that the vortex-mode vector solitons are lin- 
early unstable, we proceed with the study of the dipole 
variant. We obtain a continuous family of stationary so- 
lutions in the same parameter range. In this case, our 
tools for the linear stability analysis give no conclusive 
results, but numerical simulations of highly perturbed 
states show a periodic, stable evolution. Finally we con- 
clude and summarize the implications of this work. 

The model- We consider two incoherently interacting 
beams propagating along the direction z in a bulk, weakly 
nonlinear saturable optical medium. The model corre- 
sponds, in the isotropic approximation, to the experimen- 
tally realized solitons in photorefractive materials. The 
problem is described by the normalized, coupled equa- 
tions for the slowly varying beam envelopes, E\ and Ei. 
The normalized dynamical equations for the envelopes of 
two incoherently interacting beams can, in this case, be 
approximately written in the form 
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where Aj_ is the transverse Laplacian. Stationary so- 
lutions of Eqs. (|l|) can be found in the form E\ = 
\[&~[ u(x, y) exp(i/3iz), E 2 = y/jhv(x,y) exp(i/3 2 z), 
where Pi and 02 are two independent propagation con- 
stants. Measuring the transverse coordinates in the units 
of V/^Tj an d introducing the ratio of the propagation con- 
stants, A = (1 — /?2)/(l — Pi), from Eqs. (|T]) we derive 
a system of stationary equations for the normalized en- 
velopes u and v: 
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where /(/) = 1(1 + a/) -1 , I = u 2 + v 2 , and s = 1 - ft 
plays the role of a saturation parameter. For s = 0, this 
system describes the Kerr nonlinearity. In this paper we 
will work with intermediate values of saturation, around 
s = 0.5. 

Vortex-mode solitons.- First, we look for radi- 
ally symmetric solutions u{x,y) = u(r), w(x,y) = 
w(r) ex.p(im<j)), in which the second component carries 



a topological charge, m, and we assume that the u com- 
ponent has no charge. In this case, Eqs. (||) take the 
form: 



A T u-u + uf(I) = 0, 
A r v - (m 2 /r 2 )v - Xv + vf(I) = 0, 



(3) 



where A r = (l/r)(d/dr)(rd/dr). The fundamental, bell- 
shaped solutions with m = only exist at A = 1. In 
the remaining region of the parameter plane (s, A), the 
solutions carrying a topological charge (m = ±1) in the 
second component exist. Solutions of this type for the 
saturable nonlinearity are found and thoroughly inves- 
tigated in ||, and in Q - for the so-called threshold 
nonlinearity. The families of these radially symmetric, 
two-component vector solitons are characterized by a sin- 
gle parameter A, and at any fixed value of s, the bor- 
der of their existence domain is determined by a cut-off 
value, A c . A two-component trapped state exists only for 
A > A c . Near the cutoff point this bell-shaped state can 
be presented as a waveguide created by the u— component 
guiding a small-amplitude mode v. Away from the cut- 
off, the amplitude of the v— component grows, and the 
resulting vector soliton develops a ring-like shape. An 
example of the ring-shape vortex-mode for our model is 
presented in Figs. 2(b). An important physical charac- 
teristic of vector solitons of this type is the total power 
defined as P = P u + P v = 2tt (u 2 + v 2 )rdr, where 
the partial powers P u and P v are the integrals of motion 
for the model (|l]). The dependencies P(X), P V (X), and 
P U {X) for a fixed s completely characterize the family of 
vector solitons [a typical example is shown in Fig. 2(a)] 
for A > A c . 
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FIG. 2. (a) Power diagram for the vortex-mode vector 
soliton family for s =0.5 (solid - total power, dashed - P u , 
dotted - P v ) and (b) typical intensity profiles of the soliton 
components at A = 0.6. 

Stability analysis.- Similar to the well studied case 
of the (1+1)-D vector solitons (see, e.g., Ref. ||), the 
vortex-mode soliton is associated, in the linear limit, with 
a soliton-induced waveguide supporting a higher-order 
mode. It is therefore tempting to draw the analogy be- 
tween the higher-order (1+1)-D two-hump solitons and 
(2+l)-D ring-shape solitons. Given the established sta- 
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bility of the multi-hump one-dimensional structure in a 
saturable medium @] , this line of thought would lead us 
(erroneously!) to conclude that the vortex- mode vector 
solitons in our model should be linearly stable. 

To show that the above conclusion is wrong, we 
have performed a linear stability analysis of the two- 
dimensional vortex-mode vector solitons. Our technique 
consists in linearizing Eqs. (|l]) around the vortex so- 
lution and evolving them with completely random initial 
conditions 10 1. Usually, the solution will be a linear com- 



for large propagation distances as a rotating and radiat- 
ing pulsar state. 



bination of modes evolving with some real frequencies, /x. 
However, if the linear equation bears modes with complex 
eigenvalues \p = Re(/z) + ilm(/i)], we expect an exponen- 
tial growth of our random data, with convergence to the 
invariant space of one of these eigenvalues. This method 
allows us to extract the most unstable eigenvalue and 
its associated manifold in a similar way to the classical 
analysis of Lyapunov exponents in unstable systems. 

Our linear stability analysis has proved that, although 
saturation does have a strong stabilizing effect on the 
ring vector solitons || , all vector solitons of this type are 
linearly unstable. In Fig. 3(a), where we plot a typical de- 
pendence of the eigenvalue /i of the most unstable mode 
on the soliton parameter A for a fixed s, we see that the 
growth rate of the instability tends to zero at the cut-off 
point of the vortex mode [cf. Fig. 2(a) and 3(a)]. This 
behaviour is consistent with the inherent stability of the 
fundamental scalar soliton in a saturable medium. Else- 
where the growth rate is positive and increases when the 
increasing intensity of the vortex mode. 
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FIG. 3. (a) Eigenvalue of the leading unstable mode for 
vortex-mode solution (s = 0.5) and (b) typical evolution of 
the vortex-mode soliton near cut-off (s — 0.65, A = 0.6). 
Shown are intensity distributions of the it- and v- components 
in the (x,y) plane [—10 < (x,y) < 10]. 

In Figs. 3(a-b) we compare the linear stability analy- 
sis with dynamical simulations of the vortex-mode soliton 
near cutoff, perturbed with random noise. The instabil- 
ity, although largely suppressed by saturation, triggers 
the decay of the soliton into a dipole structure [as shown 
in Fig. 3(b)] for even a small contribution of the charged 
mode. The dipole demonstrates astonishing persistence 




FIG. 4. (a) Power diagram of the dipole-mode vector soli- 
ton family for s — 0.5 (solid - total power, dashed - P u , 
dotted - P v ), and (b) typical intensity profiles of the soliton 
components at A = 0.6. 

Dipole-mode solitons.- It is apparent from the above 
analysis that the dipole-mode vector soliton should be 
more stable than the vortex-mode soliton. We have iden- 
tified the existence domain of these solitons by solving 
numerically the stationary equations (||) to find localized 
asymmetric solutions carrying a dipole mode [as shown 
in Fig. 1(c)] in the v— component. One characteristic 
example of the dipole-mode soliton family is shown in 
Fig. 4(a) for a fixed s. Our linear stability analysis for 
these solutions does not converge to a particular value 
of lm(p,). This indicates that either the eigenvalues of 
the unstable modes, if they exist, are extremely small, or 
that the unstable modes have shapes which are only too 
weakly excited by the random perturbations. 

To obtain further information on the dynamical sta- 
bility of dipole-mode vector solitons, we have propagated 
numerically different perturbed dipole-mode solitons for 
distances up to several hundreds of z units, or diffraction 
lengths. To put these numbers into physical perspective, 
we note that, in current experiments on solitons in pho- 
torefractive materials, the typical crystal length is w 20 
mm, whereas z = 100 in our model corresponds to the 
soliton propagation length ss 20 — 40 mm. 

We have performed two types of numerical experi- 
ments. First, we have found that small perturbations 
or random noise lead to bounded oscillations of the vec- 
tor soliton, which retains its shape. This shows that the 
dipole-mode vector solitons should be stable enough for 
experimental observation. On the other hand, strong per- 
turbations, such as a large disproportion of the humps of 
the dipole or a relative displacement of the components, 
may alter the shape of the dipole. Such a perturbation, 
in the process of evolution, is typically transfered from 
one component to the other in a robust, periodic way. It 
is clear from our simulations [Fig. 5] that in this case 
the dynamics is similar to that of two beams spiraling 
together but with initial zero angular momentum Q . 
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FIG. 5. Evolution of a dipole vector soliton, (A, s) = 
(0.7, 0.5), after moving 25% of the total power from one hump 
to another. Shown are intensity distributions for the u— and 
v— components in the (x,y) plane [—7 < (x,y) < 7]. 

From the evidence above it is clear that the unstable 
modes of a dipole vector soliton, if they exist, should 
be rare and hard to excite. Calculations of the lin- 
ear stability spectrum around the exact soliton solution 
should provide a complete answer to the stability ques- 
tion, which will be the subject of our further research. 

Summing up the stability results, we can state that the 
dipole-mode vector solitons are extremely robust objects. 
According to our numerical tests, they have a typical 
lifetime of several hundred diffraction lengths and survive 
a wide range of perturbations. This is true even for vector 
states with a large power in the dipole component (P w > 
P u ), i.e. in the parameter region which cannot be safely 
reached for vortex-mode solitons, since the latter become 
unstable much sooner. 

Finally, it is important to mention that there exist 
other physical models, where stable, radially asymmet- 
ric, dipole solitary waves play an important role in the 
nonlinear wave dynamics. All these models, however, 
possess scalar, or one-component, structures. The most 
famous examples are the Larichev-Reznik soliton, a local- 
ized solution of the Charney equation for Rossby waves 
pj , and the dipole Alfven vortex soliton in an inhomoge- 
neous plasma [fl2| . Another type of dipole solitary waves 
is found as a two-soliton bound state in a non-local ]l3| or 
anisotropic nonlocal ]l4|] media due to anomalous inter- 
action of two solitons with opposite phase. Nevertheless, 
the physics behind all known dipole solitary waves and 
corresponding single-component nonlinear models differ 
drastically from the problem considered above. There- 
fore, the dipole-mode vector soliton we describe in this 
Letter is a genuinely new type of solitary waves in a ho- 
mogeneous isotropic bulk medium, a phenomenon that 
may occur in many other physical applications. 

In conclusion, we have analyzed the existence and 
stability of radially symmetric and asymmetric higher- 
order vector optical solitons in a saturable nonlinear bulk 



medium, and predicted a new type of optical soliton asso- 
ciated with the dipole mode guided by a soliton-induced 
waveguide. We have demonstrated that solitons carry- 
ing a topological charge are linearly unstable and, as a 
consequence, they may decay into dipole-mode solitons. 
There is also strong evidence of the stability of these 
dipole-mode solitons. We believe that all the effects pre- 
dicted in this Letter, including the existence and stabil- 
ity of solitons, can be easily verified in experiments with 
photorefractive media. 
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